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Abstract 

We have studied finite-sized single band Hubbard chains with Fibonacci mod- 
ulation for half filling within a mean field approximation. The ground state 
properties, together with the dc conductivity both at zero and non-zero tem- 
peratures, are calculated for such quasi-periodic Hubbard chains. While a 
reduction in the conductivity is found for strong electronic interaction or 
strong Fibonacci modulation, a competition between these two is observed 
to enhance the conductivity. The results at finite temperatures also illustrate 
some interesting features of such finite-sized systems. 

PACS No.: 71.10.Fd, 71.23.Ft, 75.30.Fd 



I. Introduction 



In recent years the study of the interplay of correlation and disorder in 
condensed matter systems has become a subject of great interest. Some re- 
cent experimental results, e.g. the metal-insulator transition (MIT) in two 
dimensional electron gas in Si-MOSFETs [1], suggest that there exists a com- 
petition between the electronic correlation and the disorder in such systems. 
Sheer complexity of the problem of considering simultaneous effects of both 
correlation and disorder attracted the attention of the researchers over the 
past few years. For example, the effect of Hubbard correlation in presence of 
random disorder has been studied by using quantum Monte Carlo techniques 
[2] and the bosonization [3]. Effect of electron-electron interaction on two- 
dimensional finite clusters of disordered spinless fermions has been studied 
by diagonalization in a truncated Hartree-Fock (HF) basis [4]. Persistent 
current in disordered mesoscopic rings with spinless fermions, threaded by a 
magnetic flux, has been calculated by using exact diagonalization technique 
for finite-sized systems [5]. Similar systems have been studied by density 
matrix renormalization group (DMRG) [6] also. It seems interesting at this 
point to study the effect of electronic interaction in a quasi-periodic lat- 
tice. This not only interpolates between the extreme cases of full grown 
order and random disorder, but also gives some insight into the electronic 
properties of quasi-crystalline superlattices which are now realizable in the 
laboratory. Since the discovery of the quasicrystals [7], the electronic prop- 
erties of these novel systems have been studied with great interest. The low 
temperature studies [8] of the electronic conductivity of the quasicrystalline 
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systems pose some interesting questions. On the other hand, recent fabri- 
cation of quasicrystalline samples with local magnetic moments [9] raise the 
issue of dealing with the electronic correlation in quasicrystalline materials. 
Consequently some works have already been carried out in this direction 
which include bosonization and weak-coupling renormalization group (RG) 
studies of interacting spinless fermions with diagonal Fibonacci modulation 
[10]. Hartree-Fock calculations for such Fibonacci-Hubbard chains have also 
been done [11] to identify the nature of the single particle spectrum and that 
of the wavefunctions in case of small electronic correlation. Recently, diago- 
nal and off-diagonal Fibonacci modulations in an otherwise periodic Hubbard 
model has been studied by a weak-coupling RG, complemented by a DMRG 
calculation [12], which makes an attempt to identify gaps in the spin/charge 
sectors. However, to the best of our knowledge, a direct calculation of the 
electronic conductivity for such a system is yet to be done. 

In view of this perspective we attempt a study of the competition between 
a quasiperiodic "disorder" (e.g. Fibonacci modulation in site potentials) and 
an on-site electronic correlation (e.g in the one band Hubbard model) in one 
dimension. The present study mainly concentrates on the calculation of the 
Drude weight [13] within a generalized Hartree-Fock approximation (GHFA) 
for studying the dc conductivity of such chains of finite lengths at half-filling. 
We have also studied the charge and spin density order parameters to char- 
acterize the ground state. It may be noted here that the study of the Drude 
weight may not be meaningful for a very large system with strong correlation 
or disorder [14]. For example, the strong on-site correlation could suppress 
the conductivity at half filling for a large chain. Therefore, our calculation is 
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restricted to small finite systems in order to resolve the issue of competition 
between correlation and disorder. In fact previous Hartree-Fock studies of 
similar systems of large sizes [11] did not attempt the calculation of conduc- 
tivity and were restricted to the limit of small correlation only (neglecting 
spin polarization) which is not expected to fully capture the essence of the 
competition between correlation and disorder. 

In section II we define the Hamiltonian and explain the calculations in 
GHFA. Section III describes the results obtained from GHFA for zero temper- 
ature at half-filling. Section IV deals with the results at finite temperature. 
Section V summarizes the present work. 

II. The Model and the GHFA calculations 

(i) Zero temperature calculation 

The Fibonacci sequence in any generation grows as: 
AB AAB AB AAB AAB AB AAB AB AAB AAB AB AAB AAB , grow- 

ing with the generating rule where A goes to AB and B goes to A in the next 
generation. For an infinite Fibonacci chain the ratio of A to B is (y/5 + l)/2, 
which is known as the golden mean. 

We take a finite chain of N sites, where N is a Fibonacci number, in the 
form of a ring that is threaded by a flux (in units of basic flux quantum = 
hc/2ne) [14]; the Hamiltonian looks like: 

H = S €iUi + t [ e% * C ia c i+la + h - c - ) + U Yl n ^ n Hi (!) 
i \ iu I i 

where q is the site energy at the i-th site; it takes on the value or e# 
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depending on whether it belongs to an A-type or a 5-type site in the Fi- 
bonacci sequence. The number operator rii a =c\ a Ci a and rii = + n^; t 
is the hopping integral between the nearest neighbour sites. The phase fac- 
tor (e 1 ^) appears because of the flux threaded by the ring. U is the on-site 
Coulomb interaction. 

On employing the GHFA the Hamiltonian is decoupled into up-spin and 
down-spin parts with the modified site energies for the up- and the down-spin 
electrons. The modified site energies are given by: 

e' ir = €i + U < > 
e-| = 6i + U < n iT > 

This decoupled Hamiltonian is then diagonalized (in a self-consistent 
manner) for both the up and the down spin parts separately. We then calcu- 
late the ground state energy E g by summing the states up to the Fermi level 
for both up and down spins corresponding to a desired filling: 

E g = J2 E ^ + J2 E nl-Uj2< n A ><rHi> ■ (2) 

i 

where E^s are the single particle energy levels obtained by diagonalizing 
the decoupled Hamiltonian for the spin species a. 

The expressions for charge density wave (CDW) and spin density wave 
(SDW) correlation functions are given by 

^) = 4lE^-^ ) K-l)(n m -l)| 

p,m 

S(Q) = hT,e iq{Rp - Rm) (n P i ~ n pl )(n m] - - n ml )\ 

p,m 

(3) 
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respectively, where q is the wave-vector of the corresponding density wave. 
R p and R m are position vectors of the p-th and the m-th sites respectively. 
The Drude weight which is a measure of the conductivity of the system is 
given by [5, 14, 15]: 

D = — — (^A (4) 

(ii) Finite temperature calculation 

Upon diagonalization of the decoupled "independent electron" Hamilto- 
nian one can calculate the partition function 

Z = Tr[e-^ H -<^\ , 

where (3 = 1/kT, k being Boltzmann's constant and T, the temperature. 
This leads to direct calculations of several thermodynamic quantities. The 
Drude weight can be evaluated from 



N ^ 



d 2 E n ^ \ i f d 2 E n i\ i 



+ lJ 



(5) 



III. Results for zero temperature 

Our results are based on calculations for different sizes of the Fibonacci 
ring e.g. iV = 34, 55, 89, 233 and 377. We started our iteration with e\ =0 
and €b = e > 0, with U > the scale of energy being chosen by t = 1. After 
achieving the self-consistent solution we find out and for each site i, 
which enables us to calculate the CDW and SDW order parameters defined 
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in (3). Eg is obtained as a function of 0, by summing over states up to the 
Fermi level for both spins. The double derivative of E(<f>) with respect to 
at 0=0 gives the Drude weight. We have carried out the calculations for 
different system sizes with a fixed value of e while varying U and vice-versa. 

It is interesting to observe how the competition between the correlation 
(U) and quasiperiodic modulation (e) decides the conductivity of the system 
(of a given size). It is also interesting to observe the dependence of this 
effect on system size. The competition between disorder and correlation is 
expected to be most prominent at half filling. This leads us to focus on this 
specific band-filling. 

Fig. 1(a) shows the variation of the Drude weight with (U/t) for different 
values of e in a half filled chain of 34 sites. In this plot one can observe that 
the Drude weight is maximum at U — for e = and the scattering due to 
the Fibonacci "disorder" (for e ^ 0) tends to lower the conductivity compared 
to the case of e = 0. However, in presence of the Fibonacci modulation in 
the site potentials (for e ^ 0), a competition between the correlation and 
the "disorder" enhances derealization and, thereby, the conductivity. As 
a result of this, there should be a gradual rise in the Drude weight with 
increasing U/t. Such a competition between the Hubbard correlation and 
the aperiodic Fibonacci modulation is indeed showing up in the plots of the 
Drude weight in Fig. 1(a) with e = 1.0 and 1.5 respectively. It can be argued 
on a heuristic ground that at U ~ 0, there will be a considerable number 
of double occupancies (two electrons with opposite spins) that would favour 
to sit on A- type sites with lower site-energies for e ^ 0; this tendency of 



6 



"pinning" would decrease the Drude weight from that of the case with e = 0. 
An increase in the value of U would tend to melt the "doublon"s and thereby 
undo the "pinning" effect. This would make the electrons more and more 
delocalized resulting in a rise of the Drude weight with increasing U/t for 
e 7^ 0. After attaining a maximum value at a certain U/t, which depends on 
the value of e (apart from the system size), the Drude weight will again fall 
with increasing U/t because of the usual suppression of the conductivity by 
the Hubbard correlation. Eventually the Drude weight becomes vanishingly 
small at a value of U/t, say (U/t) c , which also depends on e and system 
size. The general features noted above for the case of a ring of 34 sites are 
also observable in larger rings of 55, 89 sites (Fig.s 1(b) and 1(c)). Even 
for systems containing of as many as 233 or 377 sites (Fig.s 2(a) and 2(b)), 
we could easily observe similar features i.e. the qualitative dependence of 
Drude weight on e and (U/t) is independent of system size. We must make 
a comment here that the non-zero value of the Drude weight at half-filling is 
an effect of finite size. That is clearly indicated by the monotonic decrease in 
the peak value of the Drude weight with system size for a given value of e. It 
is interesting to note here that the value of (U/t) for which the Drude weight 
is a maximum decreases with increasing disorder strength (e). This feature, 
not clearly revealed for smaller system sizes can be prominently observed in 
Fig.s 2(a) and 2(b) where we find a shift of the peak in the Drude weight 
towards lower values of (U/t) with increasing e. This is so prominent in Fig. 
2(b) that for large enough e we cannot recognize the gradual rise in the Drude 
weight. The effect of disorder takes over to suppress the derealization due 
to the competition between U and e. For smaller system sizes such effect 
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will take place for larger e. On the other hand the value of (U /t) c , where the 
Drude weight drops down to zero, depends on the system size. For a given 
value of e it is observed that (U/t) c gradually decreases with the increase 
of the system size (Fig. 3). So the competitive effects of correlation and 
disorder show up in the conductivity of such systems of small sizes. The 
larger is the system size the weaker are the effects. Again for a ring of 
given length the competitive effect shows up for moderate values of e and 
U . A very high value of correlation monotonically drives the system to an 
insulating phase; scattering due to strong disorder also suppresses (even if 
U is reasonably large) the interplay of correlation and disorder (Fig. 2(b)). 
Such observations are in agreement with previous studies [12]. 

To investigate the nature of the ground state further we have plotted the 
spin and charge structure factors S(q) and C(q) (defined in (3)) respectively. 
In Fig.s 4(a) and 4(b) we have plotted S(q) as a function of q for different 
cases. In Fig. 4(a) we observe the sharp rise of the spin correlation function 
around q=it for large U (a value of U when the Drude weight is vanishingly 
small) for system sizes iV=34, 89 and 233, at e = 0. It is to be noted that 
precisely the point q=n is missing for finite odd sized systems. In the " clean" 
limit antiferromagnetic SDW phase should settle, for larger (U/t) which is 
revealed in the plot. We study the spin structure factor for e ^ in Fig. 
4(b) where we find smaller value of the q=ir peak and some wiggles at q 
7^ 7r due to an imperfect antiferromagnetic modulation caused by aperiodic 
site potentials. In presence of the Fibonacci modulation, there may arise a 
"pinning" tendency of the holes at special sites; hence the reduction in the 
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value of S(q) around q=n and the wiggles at q ^ it. The effect is more 
pronounced at higher values of e. Plots of S(q = n) as a function of (U/t) 
(Fig.s 5(a) and 5(b)) shows that the antiferromagnetic SDW grows after 
crossing a certain (U/t) that depends on the system size and on the value of 
e. These values (U/t) match well with (U/t) c values obtained in the study of 
the Drude weight. 

Fig.s 6(a), 6(b) and 6(c) show plots of the charge correlation function C(q) 
against q for different system sizes and different U's at a reasonably large 
values of e. The quasiperiodic modulation generates peaks of different heights 
at incommensurate values of q. Density waves of such incommensurate (with 
lattice spacing) wave vectors get diminished in magnitude with increasing 
U. However the position of the peaks do not change. With the increase of 
the system size the prominent peaks do not disappear but get sharper and 
sharper. Moreover additional peaks of smaller magnitude appear to indicate 
excitation of new density waves for large system sizes. This is expected 
because of hierarchical aperiodic structures of the Fibonacci generations. 

In explaining the decrease in the Drude weight with increasing "disor- 
der" we conjectured the formation of "pinned" holes and "doublons" . These 
"pinning" would follow the Fibonacci modulation and would be partly mod- 
ified by the probability of hopping between adjacent sites of comparable site 
energies. That this really takes place is evident from the plot of the density 
profile against the site index (Fig. 7) for the 34-site ring. On the other hand 
A type sites favour retaining "doublons" (with average occupancy greater 
than 1). The AA pairs, on the contrary, show a lesser preference for the 
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"doublons" to facilitate nearest neighbour hopping which lowers the energy. 

IV. Finite temperature calculation of Drude 
weight 

The Drude weight, as calculated from (5), is plotted against the tem- 
perature in Fig.s 8(a), 8(b), 9(a) and 9(b) for system sizes N = 34, 55, 89 
and 233 respectively. Cases with zero and non-zero values of e and U are 
compared in the graphs. Since the effect of correlation is expected to play an 
important role at the half-filling, the 34-site system is studied at this special 
filling (= ||) while for the 55, 89 and 233 -site chains we concentrate on 
the fillings n— 28/55, 44/89 and 116/233 respectively which are the closest 
possible cases to (i.e. one particle away from) half-filling. 

Fig. 8(a) shows the variation of the Drude weight against the temper- 
ature (kT) for system size N = 34. In absence of correlation (U = 0) the 
Drude weight starts from its maximum value (i.e. At /n) at half-filling for 
the periodic case (e = 0). There is a gradual fall in the Drude weight with 
the rise in the temperature. This is because the thermal excitations tend 
to populate the higher levels which diminish the current generated by the 
breaking of the time-reversal symmetry due to the external field. The rate 
of decrease of D, however, is too small to be clearly detected in the graph. 
However, in the presence of the "Fibonacci disorder" (e ^ 0) the scattering 
effect takes over and this makes the conductivity fall in a much sharper fash- 
ion with the increase of the temperature. On the other hand, at T = 0, 
the Hubbard correlation introduces a Mott gap [16, 17] at the middle of the 
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band below which all levels (i.e. the lower band) are occupied at half-filling. 
Therefore, the conductivity drops down to a much lower value (non-zero for a 
small finite system) in the case of U 7^ as compared to the non-interacting 
system. As the temperature increases, the excitations of the electrons take 
place across the Mott gap to populate the upper band. This results in an 
enhanced conductivity of the system. In the limit of very large temperature 
the Drude weight asymptotically approaches to that of the non-interacting 
case (e — U = 0) from below. Simultaneous presence of aperiodic "disorder" 
and the Hubbard interaction leads to a competition that reduces the Mott 
gap. This is why the value of the Drude weight for e 7^ 0, U 7^ is slightly 
higher than that in the case of e = 0, U 7^ 0. Now the Drude weight falls 
with the rise in temperature, the rate of fall being primarily dominated by 
the effect of disorder. 

More or less similar pictures evolve for other system sizes (Fig.s 8(b), 9(a) 
and 9(b)). However, a few notable differences are there. Firstly, the effect of 
the correlation becomes more pronounced with increasing system-size. This 
is why the Drude weight is vanishingly small for U/t = 1.7 (> (U/t) c ) for 
system sizes as large as N = 89 and 233 irrespective of the presence of 
disorder (i.e. for both the cases of e = and e 7^ 0). Secondly, the rate of fall 
of the Drude weight (with temperature) is enhanced for larger systems for 
which the effect of disorder-induced scattering is stronger. It is noteworthy 
that in Fig. 8(b) for N = 55 the Drude weight initially falls and then 
increases with the rise in the temperature for the ordered case with large 
correlation (e = 0, U/t — 1.7). This feature is in contradiction with the case 
depicted in Fig. 8(a) for iV = 34. This is because we have considered the 
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precisely half-filled case for N = 34 while the case of N = 55 corresponds 
to a case of one particle away from the half- filling (filling = 28/55). So, 
even in presence of the Mott gap (comparable to the finite-size gaps) the 
system behaves like a metal at low temperatures. Therefore, for the initial 
rise in temperature we observe a slow decrease in the Drude weight in Fig. 
8(b). For higher temperatures the effect of thermal excitations of electrons 
across the Mott gap dominates. This results in the ultimate rise in the Drude 
weight in this region. It is to be noted that this feature will gradually vanish 
with increasing system size for which the finite-size gaps become smaller 
and smaller in comparison with the Mott gap. Indeed, this feature is merely 
detectable in Fig. 9(a) for N = 89 while it completely disappears for N = 233 
(Fig. 9(b)). The fact that the effect of the Mott gap becomes more dominant 
for larger systems (for the same value of the Hubbard parameter U) is also 
evident from the case of e = 0.5, U = 1.7 in Fig. 9(b) i.e. for N = 233. 
Here we find a small rise in the Drude weight with temperature even in 
presence of disorder, a feature not observable in Fig.s 8(a), 8(b). For larger 
systems the introduction of a small disorder essentially drives the system to 
become an insulator which is can be sharply contrasted with that of a Mott 
insulator. In the regime of higher temperatures we find that the conductivity 
of the disorder-driven insulator (e = 0.5 U/t = 1.7) is less than that of the 
correlation-driven (Mott) insulator (e = 0, U/t = 1.7) by several orders of 
magnitude. This is, indeed, a very good illustration of differentiating the two 
types of insulators— namely, disorder-driven insulator and Mott insulator. 



V. Conclusion 
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Summarizing, we have studied finite-sized Fibonacci-Hubbard chains at 
half filling (or one particle away from half-filling) . In the present work, within 
the mean field approach of GHFA, we calculate the ground state properties 
e.g. the spin and charge correlation functions, and the dc conductivity both 
at zero and finite temperatures. We find that the dc conductivity, as obtained 
from the Drude weight [5, 14], shows a strong dependence on the competition 
between the Fibonacci modulation and the Hubbard correlation. This tends 
to enhance the conductivity at zero temperature in the regime of interme- 
diate coupling. However, in the asymptotic cases of very weak and strong 
Hubbard correlations (in presence of a reasonable modulation in the site po- 
tentials), the conductivity is reduced due to diagonal aperiodicity (driving 
towards a Fibonacci modulated charge ordering) and electronic correlation 
(bringing in antiferromagnetic fluctuations) respectively [12]. The value of 
(U/t) c , where the conductivity of the system drops down to zero, decreases 
with the system size. This suggests that the true effect of competition of 
the correlation and disorder could be appreciably observed in systems of 
smaller size (typically~100 sites). The SDW correlation strongly grows with 
Hubbard correlation at q=ir at half-filling. This indicates the dominance of 
antiferromagnetic instabilities at large U values. The charge correlation C (q) 
shows spikes characteristic of Fibonacci modulation. This charge modulation 
although reduced by the presence of finite U does not go off even at a value 
of U for which anti-ferro SDW starts dominating. There appear more and 
more characteristic spikes in C(q) as the system size increases. However, 
the spikes observed for smaller system sizes never disappear. This suggests 
that the system captures the self similar pattern of the Fibonacci hierarchy 
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and this is not masked by correlation. At finite temperatures, the Drude 
weight slowly decreases with rise in the temperature, in the periodic limit 
in absence of correlation. The presence of Hubbard correlation opens up a 
Mott gap which quenches the conductivity at low temperatures. Therefore, 
with the increase of temperature, there are excitations of electrons across 
the Mott gap to give rise to an initial increase in the conductivity. This 
semiconductor-like behaviour is in sharp contrast to the periodic case. On 
the other hand, the scattering induced by the Fibonacci disorder makes the 
picture completely different. In presence of Fibonacci modulation alone the 
Drude weight sharply decreases with temperature. This feature shows up 
even in presence of the Hubbard interaction for smaller system sizes. How- 
ever, for larger systems at very low temperatures (in presence of correlation) 
the Mott gap reduces the conductivity to such a small value that the disorder- 
driven fall in the Drude weight is almost invisible. In the high temperature 
region the behaviour of the Mott insulator and that of the disorder- driven 
insulator are in sharp contrast to each other. The results obtained for a 55- 
site system (studied for one particle away from half filling) indicate that the 
effect of variation in the band filling could give rise to some notable features. 
As revealed in the present study, the competition between electronic corre- 
lation and diagonal aperiodicity leads to interesting effects in the properties 
of such systems at both T = and T 7^ for smaller system sizes. It is, 
therefore, worthwhile to make further studies of such systems. For example, 
it would be interesting to see the effect of different types of quasiperiodic 
modulations and the effect of external fields in such systems, using similar 
mean-field approaches as well as techniques beyond the mean-field approxi- 
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mations [2, 12]. Also a detailed study of the effect of band filling away from 
the half filled case would be an interesting one at finite temperatures. 
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Figure Captions 



Fig. 1. Plot of the Drude weight D (at zero temperature) vs. U/t for different 
values of e for a Fibonacci-Hubbard ring of (a) 34 sites (with 34 particles), (b) 
55 sites (with 56 particles) and (c) 89 sites (with 88 particles). A finite value 
of D for low U/t at half-filling is a finite-size effect. For non-zero values of 
e the competition between Hubbard correlation and the diagonal Fibonacci 
modulation is clearly seen from low to intermediate values of U/t — upto a 
value (U/t) c . 

Fig. 2. Plot of the Drude weight D (at zero temperature) vs. U/t for the 
Fibonacci-Hubbard chain of (a) 233 sites (with 232 particles) and (b) 377 
sites (with 376 particles) for different values of e. The shift of the peak in 
the Drude weight towards U/t = is clearly visible for higher values of e. 

Fig. 3. Plot of (U/t) c vs. N, the system size (at/ one particle away from 
half-filling) . 

Fig. 4. Plot of the SDW correlation function S(q) (at zero temperature) 
as a function of the wave-vector q for different systems sizes (a) in absence 
of disorder (e = 0) and (b) in presence of Fibonacci modulation (e = 1.5). 
Value of the Hubbard correlation U = 3.0 (the scale of energy is given by 
t = 1.0). The peak at q = it shows antiferromagnetic spin modulation in- 
duced by large U. This correlation is reduced by the presence of Fibonacci 
"disorder" and also there appear wiggles at q ^ n. Such wiggles are merely 
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visible in (b). 



Fig. 5. Plot of S(q = n) (at zero temperature) vs. U/t for different values 
of e for the Fibonacci-Hubbard chain of (a) 34 sites, and (b) 89 sites. For 
the Hubbard model (e = 0) S(q = ir) grows with U/t at half-filling. The 
sharp rise in S(q = n) follows after a certain value of (U/t) which is in close 
agreement with the (U/t) c observed in Fig. 1. 

Fig. 6. Plot of the CDW correlation function C(q) (at zero temperature) 
vs. the wave-vector q for the Fibonacci-Hubbard chain of (a) 34 sites, (b) 
89 sites, and (c) 233 sites. There are characteristic peaks for a large value 
of the "disorder" (e = 1.5) which do not disappear with the increase in the 
correlation parameter U. 

Fig. 7. Plot of the particle density n« against the site index % for the Fibonacci- 
Hubbard chain of 34 sites at zero temperature at half-filling with e = 1.0 and 
U = 3.0 (in unit of t). The density profile shows formation of a charge den- 
sity modulation which closely resembles the Fibonacci pattern. 

Fig. 8. Plot of the Drude weight D vs. the temperature T (scaled up by the 
Boltzmann constant k) for the Fibonacci-Hubbard chain of (a) 34 sites (at 
half-filling), and (b) 55 sites (one particle away from half-filling) for different 
choices of e and U (scale of energy t = 1.0). 

Fig. 9. Plot of the Drude weight D vs. the temperature T (scaled up by the 
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Boltzmann constant k) for the Fibonacci-Hubbard chain of (a) 89 sites, and 
(b) 233 sites (for one particle away from half-filling for both the cases) for 
different choices of e and U (scale of energy t — 1.0). 
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